Superconductivity and spin correlation in organic conductors: 
a quantum Monte Carlo study 
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The d-wave pairing correlations along with spin correlation are calculated with quantum Monte 
Carlo method for the two-dimensional Hubbard model on lattice structures representing organic 
superconductors k-(BEDT-TTF)2X and (TMTSF^X. In both cases the pairing correlations for 
superconducting order parameters with nodes are found to be enhanced. The symmetry and the 
enhancement of the pairing is systematically correlated with the spin structure factor, suggesting a 
spin- fluctuation mediated pairing. We have further found that, as we deform the Fermi surface to 
make the system approach the half- filled square lattice, the coherence of the pairing saturates while 
the local pairing amplitude continues to increase. 

PACS numbers: 74.20.Mn, 74.70.Kn 



It is a great theoretical challenge to explore whether 
the superconductivity in 'exotic materials' such as the 
cuprates, heavy fermion systems, and (at least some) 
organic conductors can be encompassed into a sin- 
gle class, i.e., superconductivity mediated by spin- 
fluctuations. For more than a decade, the possibility of 
spin-fluctuation-mediated pairings has in fact been inves- 
tigated intensively both theoretically and experimentally 
for heavy fermion systems yj] and high Tc cuprates. Q 
There, one important sign of spin-fluctuation-mediated 
pairing has been the superconducting gap with nodes. 

There is now a body of accumulating experimen- 
tal evidence that organic superconductors || k-(BEDT- 
TTF) 2 X |4fU and (TMTSF) 2 X [§ also have non-s-wave 
gap like in heavy fermion or cuprate systems. The un- 
conventional pairing, along with the proximity of super- 
conductivity to antiferromagnetic or spin density wave 
state, suggest that pairing may be mediated by spin- 
fluctuation. The pseudogap-like behavior of 1/(T\T) 
above T c in k-(BEDT-TTF) 2 X is also remi- 

niscent of a similar behavior in the underdoped high Tc 
cuprates, which would suggest that electron correlation 
may play an important role there. 

Theoretically, a simplest many-body Hamiltonian to 
incorporate the electron correlation is the Hubbard 
model. Some analytical calculations have supported 
spin-fluctuation mediated pairing in the Hubbard model 
on lattices representing k-(BEDT-TTF) 2 X or 
(TMTSF) 2 X ]l7j . However, numerical evidences sup- 
porting such a possibility have yet to come. 

Thus the purpose of the present paper is to ex- 
plore numerically the pairing correlation in the Hubbard 
model for lattices representing k-(BEDT-TTF) 2 X and 
(TMTSF) 2 X, with a special attention payed to whether 
the pairing is linked with the behavior of the spin corre- 
lation. Comparing the results for the two cases, we find 
that the symmetry of the pairing is indeed determined 
by the the peak position of the spin structure factor, 



FIG. 1. The hopping parameters considered. 



which supports spin-fluctuation mediated pairing. We 
further deform the Fermi surface from the case corre- 
sponding to k-(BEDT-TTF) 2 X to find that the coherence 
of pairs saturates as the system approaches the half-filled 
square lattice (i.e., where, at least for large enough U, 
the system becomes a Mott insulator with antiferromag- 
netic order). Here we adopt the ground-state, canonical- 
ensemble quantum Monte Carlo (QMC) method. 

We start from the Hubbard model on a two- 
dimensional lattice (Fig j^) , 

^ = — y , [^( c -E,j/,(T C a;+l,J/,o") ~)~ tyi c x,y,o- c x,y+l,cr) 



+ t'/( c l,y,a c x+l,y+l,a) + ^(4,^^-1^+1^) + h.C 



(1) 



Here, (x, y) is the coordinate of a site with the lattice con- 
stant taken to be unity, and periodic boundary condition 
is assumed. 

While QMC method has been widely used to investi- 
gate the Hubbard model on the 2D square lattice, en- 
hanced pairing correlation had eluded detection. Re- 
cently, however, the present authors identified this as 
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coming form the fact that the pair-scattering processes 
that produce superconductivity has a very small energy 
scale of O(0. Oli) or less [jl8lfL9|: if one takes a closed shell 
condition (a parameter set with no ground state degen- 
eracy for U = 0), the energy gap between the highest 
occupied levels (HOLs) and the lowest unoccupied levels 
(LULs) for U — is as large as ~ 0(0. It) for tractable 
system sizes, so that the effect of the low energy pair 
scatterings would be smeared out. On the other hand, 
QMC calculations with an open-shell condition, in which 
the effect of low-energy pair scatterings is expected to be 
incorporated, suffer from numerical difficulties such as 
the negative-sign problem. An open-shell configuration 
is ill-conditioned also in the sense that there are a finite 
number of levels within an infinitesimal distance from the 
Fermi energy, so that one can suspect an enhancement of 
pairing correlations, if any, may be due to such an effect. 

We have circumvented the problem by taking slightly 
different values of t x and t y to lift the degeneracy between 
wave numbers (fcijfe) and (k2,k±), and put the Fermi 
level in between those levels. This way we can prevent the 
low-energy pair scattering processes from being masked 
and at the same time take a closed shell configuration. If 
we take typically t y /t x = 0.999, which gives HOL-LUL 
gap A of < O.Oli, and U/t — 1, we can achieve conver- 
gence with respect to the projection imaginary time r in 
the QMC without running into serious sign problem, pcj ] 

We calculate the pairing correlation functions, 
P(r)= Yl (0^(x + Ax,y + Ay)0(x,y) 

\Ax\ + \Ay\=r 

+0(x + Ax,y + Ay)0\x,y)) (2) 
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FIG. 2. QMC result for the d x 2_ y 2 pairing correlation as a 
function of real space distance for 134 electrons in 12 x 12 sites 
(n = 0.93) with t x = 1, t y = 0.999, t[ = 0.70, t) = -0.11, 
U = 1 (solid circles) and U — (dashed line). The inset 
shows the HOLs and LULs within 0.01 to the Fermi energy. 
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FIG. 3. QMC result for the spin structure factor as a 
function of the wave vector. The parameters are the same as 
in Fig|. 



8 X ,5 y ,cr 

(3) 

which includes the conventional d x 2_ y 2 symmetry (with 
5 X = 5 y = 1, corresponding to the order parameter 
proportional to /(k) = cos(k x ) — cos(k y ) in fc-space), 
but written here in a general form. Hereafter we define 
r = \Ax\ + | Ay | as the real space distance. To give an 
above-mentioned insight, we also calculate the spin struc- 
ture factor, 



We first look into the case where the Fermi surface, 
represented by HOLs and LULs, resembles those of n- 
(BEDT-TTF) 2 X materials. It has been proposed that k- 
(BEDT-TTF) 2 X can be modeled by a single band Hamil- 
tonian on an anisotropic triangular lattice with t x = t y 
^ ti > t', and n = 1 (half-filled). This is be- 

cause (i) the dimcrization of the BEDT-TTF molecules 
is so strong that we can consider a dimer as a single site 
as far as the low-energy excitations are concerned, and 
(ii) the hopping parameter in the c direction alternates 
between two slightly different values, but the difference 
is small enough to be neglected as a first approximation. 

First we take 134 electrons in a 12 x 12 (band filling 
n = 0.93) lattice with t x = 1, t y = 0.999, ti = 0.70, 
t' f = -0.11. We have taken t y = 0.999 for the slight lift 
of degeneracy mentioned above. The reason for taking 
small but negative t', and the slight deviation from n = 1 
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(to which the k-(BEDT-TTF) 2 X system corresponds) is 
to distribute HOLs and LULs on the Fermi surface as 
uniformly as possible. In FigJ|, we show the result for 
the d x 2_ y 2 pairing correlation P x 2_ y i(r) as a function of 
the real space distance r for U = 1. It can be seen that 
the pairing correlation for U — 1 is enhanced over that 
for C/ = 0, especially at large distances. 

If we look at the spin structure factor in fc-space in 
Fig.||, a broad peak around (n, n) is seen. Such an an- 
tiferromagnetic spin-fluctuation enhances repulsive (in 
momentum-space) pair scattering processes with momen- 
tum transfer (it, ir). Repulsive pair scattering processes 
from ~ ±(0, it) to ~ ±(tt, 0) (and vice versa) favor a su- 
perconducting order parameter that takes the maximum 
of its absolute value with different signs around these 
two areas. Thus, the enhancement of the d x i_ y i pair- 
ing correlation along with the peak of S*(q) at (it, it) are 
consistent with the spin-fluctuation mediated pairing. 

We have further investigated the link between d x 2_ y 2 
pairing correlation and the spin correlation systemati- 
cally by deforming the Fermi surface. In Fig.||, the en- 
hancement of the pairing correlation is plotted against 
S(Q) with Q = (7r,7r) for three cases I, II, and III. The 
main difference among the three cases is in the value of f ( 
as given in the caption, while other slight differences are 
due to technical reasons. Namely, we are going from the 
isotropic triangular lattice (I) to a case close the square 
lattice (III) via case (II) (corresponding to Figj^) for n 
fixed around unity. S(Q) increases as we approach the 
half-filled square lattice. In Fig^, the pairing is probed 
from two quantities: the pairing correlation (measured 
from the noninteracting case) summed over large dis- 
tances, plr = Er>4l[- P 'i( , ')lf=i ~ [Pd(r)]u=o}, and the 
local amplitude of pairing, p = [P d (0)]u=i ~ [Pd(0)]u=o- 

Fig.^ shows that po grows hand in hand with S'(Q), 
endorsing the spin-fluctuation-mediated formation of 
d x 2_ y 2 pairs. On the other hand, plr, the enhancement 
of the long-range part of Pd(r), which measures the co- 
herence of the pairs, also grows with 5(Q) as we go from 
(I) to (II), but it saturates between (II) and (III), sug- 
gesting that the coherence does not necessarily grow with 
the local pair amplitude. 

To show that such a correlation between the d-wave 
pairing and the spin correlations is not accidental for 
the cases studied above, we move on to the square lat- 



tice (t/ 
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0) for various values of the band fill- 



ing n. In Fig|| we show a plot similar to Fig.[| for the 
square lattice, along with the plot against n in the in- 
set. For intermediate densities, plr, Po, and 5(Q) all 
grow as n approaches the half- filling (n = 1). If we come 
too close to half-filling, however, plr becomes saturated, 
while S'(Q) and po keeps growing. This is even more 
clearly seen at exactly half-filling, where plr is strongly 
suppressed (inset of Fig.||), while po continues to grow. 
The result that, when S(Q) becomes too large, the pairs 
are formed but their coherence (jplr) stops to grow un- 



like their local amplitude po might have some relevance 
to the normal-state pseudogap behavior observed close 
to the superconducting-antiferromagnetic boundary in k- 
(BEDT-TTF) 2 X as well as in the underdoped high T c 
cuprates. 



0.20 



0.10 




0.08 



^3 



0.04 



0.00 

0.10 0.20 0.30 

S(Q) 

FIG. 4. A plot for plr and po 
Q = ("",1") for parameter sets (I) n ■ 
1.00, t) = 0, (II) n = 0.93, t 



0.00 



-0.11 (corresponding to Figji]), (III) 



0.40 



against S(Q) with 
= 0.96 t y = 1.001, 
0.999, t{ = 0.70, 
0.94, 



t v = 0.999, t{ = 0.20, t) = 0, all with 12 x 12 sites, t x 
and [7 = 1. 
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FIG. 5. A plot similar to Fig.^| for a 10 x 10 square 
lattice with (I) 34 (n = 0.34) (II) 66 (n = 0.66) (III) 92 
(n = 0.92) electrons. The hopping parameters are fixed at 
t x = 1, t y — 0.999, t 1 J = t'\ — 0. The inset shows the band 
filling dependence of po and Plr, including the results for 
half-filling (n = 1). 
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In order to confirm spin-fluctuation mediated pair- 
ing to a wider extent, we next explore the case where 
the spin structure factor is not peaked around (it, it), 
by studying the case with a Fermi surface representing 
(TMTSF)2X. If we neglect the weak dimerization along 
the stacking direction, these materials may be modeled 
by a single band Hamiltonian on a strongly anisotropic 
two-dimensional lattice with t x /ty ~ O(0. 1) (i^ = t', = 0) 
and n — 0.5 (quarter-filled). 0,^4) Here we take 72 par- 
ticles in 12 X 12 (n = 0.5), t x = 1, and t y = 0.212, for 
which the Fermi surface is as depicted in the lower inset 
of FigJ|(a). 

For this Fermi surface consisting of warped parallel 
lines around k x ± tt/4, the nesting vector is (tt/2,tt), so 
that the spin structure factor (lower inset of Fig.g(b)) 
has a peak there. J2{| Then, pair scattering processes 
from ~ ±(0, 7r) to ~ ±(ir/2, 0) (and vice versa) would be 
enhanced, which favors a superconducting order param- 
eter proportional to /(k) = cos(2fc x ) — cos(k y ) (5 X = 2, 
S y = 1 in Eq.(||)) rather than the ordinary d x 2_ y 2 pairing 
with /(k) = cos(k x ) — cos(k y ). In fact, the possibility of 
such a pairing was pointed out in an RPA calculation. 



In Fig.^J, we can see that the pairing correlation for 
/(k) = cos(2fc x ) — cos(ky) is indeed enhanced, while 
d x 2_ y 2 pairing correlation is not. The result provides 
another indication that the pairing symmetry is dictated 
by the dominant spin-fluctuation. 
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FIG. 6. A plot similar to Fig.| for 72 electrons in 12 x 12 
sites (n = 0.5) with t x = 1, t y = 0.212, t{ = t'/ = 0. Pairing 
correlations for /(k) = cosk x — cosk y (a) (ordinary d x 2_ y 2), 
and for /(k) = cos 2^ — cosk y (b). The upper inset in each 
figure schematically depicts the pairing (•: — , o: +) in real 
space. The lower inset in (a) shows HOLs and LULs within 
0.01 to the Fermi energy, while the spin structure factor is 
displayed in (b). 



To summarize, the present result supports the spin- 
fluctuation mediated anisotropic pairing scenario in or- 
ganic superconductors, by the same token as in the high 
To cuprates discussed in pq , pj] . The enhancement of 
the local pair amplitude grows with S'(Q), while the off- 
diagonal long-range order (plr) stops to do so as the 
system approaches half-filled square lattice. It is an inter- 
esting future problem whether this is caused by the mag- 
netic ordering or the metal-insulator transition, which is 
difficult to discern at present, since both of them occur 
at half-filled square lattice (at least for sufficiently large 
U). The possibility of the triplet pairing proposed in ref. 
fll4f also poses an intriguing issue. 

Numerical calculations were performed at the Super- 
computer Center, Institute for Solid State Physics, Uni- 
versity of Tokyo, and at the Computer Center of the 
University of Tokyo. K.K. acknowledges support by the 
Grant-in- Aid for Scientific Research from the Ministry of 
Education of Japan. 
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